Cross-Molecular Coupling in Combined Photoassociation and Feshbach Resonances 
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We model combined photoassociation and Feshbach resonances in a Bose-Einstein condensate. 
When the magnetic field is far-off resonance, cross coupling between the two target molecules- 
enabled by the shared dissociation continuum-leads to an anomalous dispersive shift in the position 
of laser resonance, as well as unprecedented elimination and enhancement of resonant photoassoci- 
ation via quantum interference. For off-resonant lasers, a dispersive shift and quantum interference 
appear similarly in resonant three-body Feshbach losses, except that the Feshbach node is tunable 
with intensity. 
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/nirorfwcfion. -Photoassociation occurs when a pair of 
atoms absorb a laser photon and thereby jump from the 
free-atom continuum to a bound molecular state [l|. At 
the turn of the last century, it was predicted that pho- 
toassociation could convert a condensate of atoms into 
a condensate of molecules [2] which, in turn, raised the 
question of a rate limit on atom-molecule conversion in 
a condensate. The rate limit on photoassociative atom- 
molecule conversion arises either from two-body unitar- 
ity [3| , or many-body rogue photodissociation to noncon- 
densate atom pairs [J, l5|. In the unitary limit, the De- 
Broglie wavelength sets the length scale i = Kd , whereas 
in the rogue limit the length scale is set by the interpar- 
ticle distance (. = p~^^^, where p is the system density. 
In either case, the fastest a molecular condensate can be 
created is ~ ml? /h. 

Early condensate photoassociation experiments fo- 
cused on bulk molecule formation [6], but next- 
generation experiments turned to the strongly interact- 
ing regime and the rate limit on atom-molecule for- 
mation [Tl, Isj. Experiments with a Na condensate at 
NIST were thwarted by strong dipole forces [7|, limit- 
ing the available laser intensity. Despite an intensity 
^ 1 kW/cm^, the rate limit remained out of reach, 
whereas an intensity-dependent redshift of the photoas- 
sociation resonance was measured to be consistent with 
previous theory [9| and nondegenerate experiments [lO|. 
The experiments at Rice focused on ''Li [8|, and a laser 
intensity ~ 80 W/cm^ was sufficient to achieve a rate 
limit consistent with unitarity. However, the system was 
only borderline quantum degenerate, and the rogue limit 
could not be ruled out. 

To probe the rate limit deeper, a Feshbach resonance 
was combined with photoassociation [11]. Also known as 
magnetoassociation, a Feshbach resonance occurs when 
one atom from a colliding pair spin flips in the pres- 
ence of a properly tuned magnetic field JI4I and, sim- 
ilar to photoassociation, the pair jumps from the free- 
atom continuum to a bound molecular state. A mag- 



netoassociation resonance enables a tunable interatomic 
scattering length [12, I3I, Il4l . and can alleviate conden- 
sate instability problems [l5|. Moreover, the Feshbach 
resonance is known to enhance laiilll or suppress 17 1 
photoassociation losses, and enhancement in particular 
could discern between the rate limits. In addition to 
an anomalous shift of laser resonance that is blue (red) 
for magnetic fields below (above) resonance, the latest 
experiments on Feshbach-assisted photoassociation of a 
condensate observe a rate constant that essentially van- 
ishes (~ 10~^^ cm'^/s) at a particular below-resonance 
magnetic field and saturates on-resonance at an unprece- 
dented value (-- 10^"^ cm^/s) [11]. 

This Letter develops a simple analytical model of com- 
bined photoassociation and Feshbach resonances. Models 
presently are based on the idea of photoassociation with 



a Feshbach-tunable interatomic scattering length [11 
id Il7| , and apply only for off-resonant magnetic fields 



Herein the basic quasi-continuum model is akin to two- 
color laser spectroscopy, valid in principle for simulta- 
neous resonance. Moreover, while our analytical model 
ultimately encompasses models of photoassociation with 
a Feshbach-tunable scattering length [111 . Il6l . Il7| . and 
agrees reasonably with observation [llj, we also foretell 
results for magnetoassociation with a laser-tunable scat- 
tering length. An experimental distinction between the 
unitary and rogue limits therefore remains elusive. 

Morfe/.-Consider N atoms that have Bose condensed 
into, say, the zero momentum plane-wave (fik = 0) state 
|0). Photoassociation and the Feshbach resonance then 
couple atoms in the state |0) to diatomic molecules of 
zero momentum in the states |1) and |2), respectively. As 
per Fig.flja), this is the F-system familiar from few- level 
laser spectroscopy. Annihilation of an atom (molecule) 
of mass m [M = 2m) from the atomic (ith molecular) 
condensate is represented by the second-quantized op- 
erator oq = a {hi). This theory is the simplest, and 
molecules dissociate only back to the level |0). To be 
more complete, molecular dissociation to noncondensate 
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FIG. 1: (color online) Few-level scheme for a condensate 
tuned nearby a combined photoassociation and Feshbach res- 
onance, (a) Basic three-level scheme, where a photoassoci- 
ation and Feshbach resonance couple the atomic condensate 
|0) and molecular the condensates |1) and |2), respectively, 
(b) A quasicontinuum accounts for dissociation to nonconden- 
sate levels, (c) Eliminating the noncondensate levels leads to 
an effective V^-system, where the virtual continuum couples 
the two molecular states, and where the detunings Si include 
the free-bound redshift. (d) When the system is far from one 
resonance, magnetic or laser, the off-resonant molecular state 
can also be eliminated, leaving an effective two-level system, 
where the detuning u includes an anomalous Stark-shift. 



levels should also be included [Fig. [Ub)]. These lev- 
els must be considered because a condensate molecule 
need not dissociate back to the atomic condensate, but 
may just as well create a pair of atoms with equal-and- 
opposite momentum, since total momentum is conserved. 
So-called rogue [J, l5|i or unwanted 18], dissociation to 
noncondensate modes therefore introduces the operators 
a±k, as well as the kinetic energy fek = h^k"^ /2vir of an 
atom pair of reduced mass rrir = m/2. 

To obtain mean-field equations, the Heisenberg equa- 
tion for a given operator is derived from the Hamiltonian 
(not shown), ihx = [x,H], and all operators are subse- 
quently declared c-numbers. In a minimalist model, x 
represents either the atomic (ith molecular) operator Ok 
{hi), or the anomalous density operator Ak = akfl-k, 
where A\^ arises from rogue dissociation to nonconden- 
sate atom pairs of equal-and-opposite momentum. Con- 
verting summations over k to integrals over frequency e 
introduces the characteristic frequency LOp ~ hp^'^ /2'mr. 



The corresponding mean-field model is (i = 1, 2) 



-^ia*bi — 5120*^2 



(la) 



ih. 



~S'A 



^n,a^ 



kJdsV^Me)A{e), (lb) 



iA{e) = eAie) ~ n, Me)b, ^ ^2 f2{e)b2. (Ic) 

The amplitudes are of unit order, the ith atom-molecule 



coupling is fi^ ex y/p, the rogue dissociation coupling is 

£,i = ^{/{in'^ujp ), the frequency dependence of the ith 
atom- molecule coupling is fi{s), and the tunable binding 
energy of the ith molecular state is hS^ . Lastly, sponta- 
neous decay of the photoassociation molecule has been 
included as 5[ = S[ — iTs/2, and spontaneous decay of 
the Feshbach molecule [19| has been neglected. 

It remains to model the continuum by specifying the 
shapes fi{e). If a condensed-matter-type universality is 
to be the driving paradigm then the only lengthscale 
in the strongly-interacting problem is the interatomic 
distance, and a single theta- function, Ji — Q(e — j3i) 
with j3i = ojp, could be employed. However, universal- 
ity has already failed in magnetoassociation of ^Li [20|, 
where the proper length scale was the size of the Fesh- 
bach molecule, and we therefore allow for different length 
scales with the Lorentzian ff — 1/(1 + 4e^//3|). 

Effective Few-Level Systems. -Here our focus will be on 
the situation where only one field, laser or magnetic, is 
resonant, and the other is off resonant. Numerically, 
Eqs. H]) are stiff for large detunings ((5^ 3> ^i), and 
an analytical solution is therefore enabled by eliminat- 
ing the rogue amplitude using A « 0; Eq. (fTc| yields 
A{e) = {flifibi + fl2f2b2)/£, which is put into Eq. (fTb| . 
and dissociation is then re-introduced phenomenologi- 
cally. This simple approximation is related to the Fermi 
Golden Rule and, like the Golden Rule, arguably reveals 
the essential physics. The model ([T]) becomes 
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-flia*bi — Q,2a*b2, 
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where Si = Si - ^(F, + Fi)/2, S2 = S2 - ^^2/2, 
and where fl^ ^ftifh/iin^cjl^^) J defiie)f2ie)/^ « 
{r|illfl2/4:nlJp)^y(3l/uJ^; note, 77 — 0.8346 is leftover 
from a hypergeometric function in the limit of point- 
like Feshbach molecules (/?2 ^ /3i). The detunings 
have the usual [9| redshift Si — S'^ — E^, where S^ = 
nl/{8TrW/^)Jdeffie)/^ = n^/ilGTTiVp),/^^;. The 
dissociation rates are F^ = fl^ /(STrLUg) y/ej /lOq, with hei 
the energy of the atom pair [2l| . Incidentally, Ei and F^ 
are the real and imaginary parts of the Lamb shift from 
laser spectroscopy. Last, the shared continuum acts like 
a virtual state that couples the photoassociation and Fes- 
hbach molecules with Rabi frequency fia [Fig. [IJc)]. 

Now, starting with an off-resonant magnetic field 
((52 > r22,r2), we take 62/^2 « in Eq. (gc]); hence, 
&2 = {^20,^ + ^3bi)/2S2 is substituted into Eqs. ^, 
yielding the two-mode system [Fig. [TJd) with |i) = |1)]: 
ia — — xa*&i, ibi — {v — iV /2)bi — jxa^, where x = 
Ui + 2/:^202/^3, ^1 - ^1 - -C^a, and F = F, -H Fi -H CT2, 
where C = ri§/4|^2|^. Note that, although the usual 
singularities are absent for S2 — 0, the model is still 
only valid far-off resonance {S2 ^ fi2,r2). Beyond the 
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usual Feshbach mean- field shift, \a 
has been neglected compared to F 
molecular" shift of the laser resonance 



ri|(52/2|52| , which lated to the magnetic field by M2 = {B 



ll| , we find a "cross- 
S3 = —Cd2, that 
is blue (red) below (above) the Feshbach resonance. Just 
as surprising is the fact that effective photoassociation 
ceases, i.e., x = 0, for 62 ~ —512^3/2^21. Borrow- 
ing again from laser spectroscopy, this node arises from 
destructive interference between direct photoassociation 
and photoassociation via the Feshbach molecular state. 
Similarly, constructive interference occurs above reso- 
nance, but observation is complicated by condensate in- 
stability due to a negative resonant scattering length 11 1 . 

On the other hand, we may also consider an off- 
resonant laser {61 ^ rii,rs,Fi), which leads to essen- 
tially the same two-mode system [Fig. [Ijd) with |i) = 
|2)]: ia = -xa*&2, ih = (^^ - «r/2)62 - 5X0^, ex- 
cept that X = ^2 + ^LSiili/il^, V — 82 — C5i, and 
F = Fs 4^ Fi + CT2, where C = f7i/4|^ip. In addi- 
tion to a dispersive shift C5i, we find that x = and 
atom losses due to magnetoassociation, i.e., three-body 
losses near a Feshbach resonance, vanish for laser detun- 
ing (5i « — J7iri3/2r22. The Feshbach node of course arises 
from destructive interference between direct magnetoas- 
sociation and magnetoassociation via the photoassocia- 
tion state. Moreover, whereas the strength of a given Fes- 
hbach resonance-and thus the detuning position of the 
above photoassociation node-is fixed, the laser intensity 
can be varied. For arbitrary red detuning —61 > 10 F^ 
and f2x(3) — J7i(3)V/, the intensity position of the Fesh- 
bach node is /„ = —26ift2/ {^1^3)- 

Comparison to Observation.~To compare with the lat- 
est observations [ll[, we first need the rate equation for 
losses from the atomic condensate l4|: Pq = —pKP^, 
where Pq = \a\^ and pK = ^X^T/li^"^ + T^/A) defines 
the rate constant K. Next, the photoassociation param- 
eters ill, ei, and pi are approximated by comparing to 
previous experiments [8|, where the natural molecular 
linewidth is Fs = 12 x 27r MHz, using the definitions 
Qi = fli^/l, ei = pl/2mr and /3i = h/{2mrL^). In 
the low-intensity fimit, we find fliy/7o = 290 x 27r kHz, 
with the saturation intensity Iq — 28 W/cm^ defined 
by Fi/Fs = I/Iq. From the high-intensity limit, the 
characteristic momentum is pi = 2.21 h/X, roughly twice 
the photon recoil momentum for light of wavelength 
2ttX — 671 nm. From the lightshift, Lf3 = 116ao is 
comparable to the classical size of the photoassociation 
molecule [8|, [ill (oq is the Bohr radius). The main pa- 
rameters for the Feshbach coupling are the product of the 
zero-field scattering length and the resonance width ll| , 
jabklAs = 1.6 nm-G, and the difference in magnetic 
moments between the Feshbach molecule and the free- 
atom pair |ll|, A^ = 2/io (the Bohr magneton is /iq), 
so that Q2 — ^87r/9|abk|ABA^/m = 127 x 2-k kHz. 
The cross- molecular coupling is then 5^3 = fJav/, where 
^^3^ = 138 X 27r MHz. The Feshbach detuning is re- 



Bo )A^, where 



Bq = 736 G [ll| locates resonance. The sole unfixed 



parameter is the kinetic energy of the magnetodissoci- 
ated pair, hi.2 = p\l2mT. For a zero-temperature ho- 
mogeneous system of point-like Feshbach molecules, the 
reasonable ansatz for the length scale is p2 = hp^l^ . 

Results are shown in Fig. [5] for condensate density 
p = 10^^ cm~^ and a laser tuned to lightshifted reso- 
nance. The net lightshift per unit intensity is E' = 
(S1+S3)// = (F,/2Jo)\/^^+£<52 with£ = ni/4|52h 
for B = 732 G, the net hghtshift is blue S' = -13 x 
27r MHz/YW/cm^), in reasonable agreement with obser- 
vation |11|. Far from the Feshbach resonance, laser- 
resonant losses approach those for photoassociation alone 
[panel (a), solid line]; but, as the Feshbach resonance is 
approached, the saturation intensity decreases and the 
rate limit increases [panel (a), dashed and dash-dotted 
lines]. From the definition F = Fg + Fi + /IF2, satu- 
ration sets in when Fi + CT2 ^ F^, which translates 
into a saturation intensity /1//0 = 1/(1 + 'C/oF2/Fs); 
for B = 732 G, we find h = 12 W/cm^. As for the rate 
limit, it is roughly the rate for converting atoms into Fes- 
hbach molecules, Rq ~ 2ri2/F2 = IGnuip. This estimate 
is best near resonance: for B = 732 G and 735.5 G, ex- 
act results from the definition of K are Ro/16Trujp — 0.35 
and 0.84, respectively. The B — 732 G results are in 
reasonable agreement with Ref. [l]|. Lastly, losses cease 
for S2 ~ -fl2^3/2rJi, or B = 714 G [Fig. H^b)], which 
is intensity independent but depends on the classical size 
of the photoassociated molecule through flsoc L/3, again 
in reasonable agreement with experiments [lll |. 

For a magnetodissociation momentum set by the inter- 
particle distance, theory compares reasonably with obser- 
vation. However, since the off-resonant size of the Fesh- 
bach molecule is roughly the resonant scattering length. 
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FIG. 2: Laser-resonant rate constant vs. (a) intensity and 
(b) magnetic field for photoassociation of a ^Li condensate 
tuned nearby a Feshbacii resonance. In panel (a), the rate 
constant saturates at a limit and intensity that depends on 
magnetic field, where the solid (dashed, dash-dotted) line is 
for B = 728 G (732 G, 735.5 G). In panel (b), the rate con- 
stant for 7 = 10 W/cm^ approaches the result for photoasso- 
ciation alone at low fields, essentially vanishes for B — 714 G, 
and rises to unprecedented values near resonance. Note that 
the model is on the edge of validity at B = 735.5 G. 



and since p~^l'^ ~ arcs(732 G) [11|, the correct length 
scale is somewhat ambiguous. We expect that mod- 
elling near-Feshbach resonance experiments will require 
the classical size of the Feshbach molecule. Finally, in 
that the saturation limit is set by the rate for magnetoas- 
sociation alone, Vo-Khp^l^ jvn^ these results agree with the 
rogue model [5| up to a dimensionless constant. However, 
the observed unitary limit 1]| coincides with the rogue 
model herein at -B = 732 G, and a definitive distinction 
remains elusive. Also, future work is needed to deter- 
mine if the rate of Feshbach-assisted photoassociation of 
a condensate saturates or maximizes. 

Before clo sing , we make the connection to existing 
models ll|, ll7[ of photoassociation near a Feshbach 
resonance. For an off-resonant magnetic field (recall, 
(52 ^ ^2,r2), the Feshbach-resonant interatomic scat- 
tering length is defined 47r?iparcs/"i — ~Vi\l2bi, so that 
the effective coupling becomes % = r2i(f — rja^-cs / L p) . 
The photoassociation node occurs when the Feshbach- 
resonant scattering length equals the classical size of the 
photoassociated molecule, Orcs ~ Lp. Not incidentally, 
what appears to be the correct definition of ares leads, 
ultimately, to the IBtt in the rate limit. Similarly, the 
net lightshift can be written S = Si(l — 2rj^ a^cs / L p) ^ 
which crosses zero, i.e., blue to-and-from red, at about 
the same detuning position as the photoassociation node. 
Of course, the atom-atom scattering length can also be 
tuned with off-resonant photoassociation [J, [oj, \2f^ , and 
analogous results apply to the Feshbach nodes. 

Conclusion.-We have reported a general model for a 
Bose-Einstein condensate near a combined photoassocia- 
tion and Feshbach resonance. When the magnetic field is 
far-off resonance, cross molecular coupling between the 
two target molecules leads to an anomalous dispersive 
shift in the position of laser resonance. Moreover, even 
though the magnetic field is far from resonance, and the 
probablility of forming Feshbach molecules is miniscule, 
this cross-coupling can eliminate or enhance resonant 
photoassociation via quantum interference. Unfortu- 
nately, a definitive experimental distinction between the 
unitary and rogue limits for photoassociation of a con- 
densate remains elusive. Nevertheless, for a far-detuned 
laser, a similar shift and interference arises in resonant 
Feshbach losses, except that the interference is tunable 
with laser intensity. These results are typical of two-color 
laser spectroscopy, despite the presence of only a single 
laser. 
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